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The motion of micro- and nanomechanical resonators can be coupled to electromagnetic fields. 
Such optomechanical setups allow one to explore the interaction of light and matter in a new regime 
at the boundary between quantum and classical physics. We propose an approach to investigate 
non-equilibrium photon dynamics driven by mechanical motion in a recently developed setup with 
a membrane between two mirrors, where photons can be shuttled between the two halves of the 
cavity. For modest driving strength we predict the possibility to observe an Autler-Townes splitting 
indicative of Rabi dynamics. For large drive, we show that this system displays Landau-Zener- 
Stueckelberg dynamics originally known from atomic two- state systems. 
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Landau-Zener (LZ) transitions [U |2] are essential to 
the dynamics of many physical systems. In the usual 
model, a parameter in a two-state Hamiltonian is swept 
through an avoided level crossing where the two bare 
eigenstates |1) and |2) hybridize. When the parameter is 
changed at a finite speed, the system may undergo a LZ 
transition into the other eigenstate. Beyond this stan- 
dard LZ problem, the dynamics becomes more elaborate 
if repeated transitions are take into account. For a peri- 
odic modulation of the parameter, the first LZ transition 
splits the state into a coherent superposition +/3|2). 
Due to the difference in energy, the system afterwards 
accumulates a relative phase between states |1) and |2). 
Thus, when returned to the avoided crossing, the system 
undergoes quantum interference with itself during the 
second LZ transition. This leads to interference patterns 
for the state population, so called Stueckelberg oscilla- 
tions [3]. Originally, Landau-Zener-Stueckelberg (LZS) 
dynamics was studied in atomic systems [H |5] |6]. Re- 
cently, the concept has been applied to superconduct- 
ing qubits [7 . Currently, there is growing interest in LZ 
and LZS dynamics concerning topics such as state prepa- 
ration and entanglement [HI |9], cooling or qubit spec- 
troscopy [10 . Another rapidly evolving area of research 
is optomechanics (see [llj for a recent review and further 
references). Optomechanical systems couple mechanical 
degrees of freedom to radiation fields. This provides new 
means to manipulate both the light field and the mechan- 
ical motion. Apart from the hope to eventually explore 
the quantum regime of mechanical motion, there have 
been several studies of the complex nonlinear dynamics 
of these systems [H [El [H |15] . 

Here, we propose an approach to observe dynamics of 
the light field in a driven optomechanical system, in the 
form of LZS oscillations. We note that there exist some 
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Figure 1: (a) Setup: a dielectric membrane couples two modes 
aL, clr inside a cavity. The left hand side is excited by a laser 
(jJl while the transmission to the right is recorded, (b) Op- 
tical resonance frequency as function of displacement: the 
membrane's displacement linearly changes the bare mode fre- 
quencies (dashed). Due to the coupling there is an avoided 
crossing of the eigenfrequencies uj± (black). The membrane 
is driven, with x{t) = Acos(Qt) + xq (blue). 

purely optical setups [16j [17] that have mimicked quan- 
tum two-state and standard LZ dynamics (but not LZS 
oscillations). In the optomechanical setup analyzed here, 
the mechanical motion of a membrane placed between 
two fixed mirrors is driven such that the resulting motion 
shuttles photons between the two halves of the cavity. A 
setup of this kind was recently realized in [181 IE] • 

We consider two cavity modes coupled by a dielectric 
membrane placed in the middle between two high-finesse 
mirrors, see Fig. [l^. The system Hamiltonian reads 

^ f. . ^ ^ / x{t)\ . ^ 
Hsys = nwo II p 1 a^aL + hcjo I H p 1 ct^RCtR 

+ hg (d[dR + d^RdLj + ^drive + ^decay. (1) 
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Oj^aL and a\^aR are the photon numbers for the two 
optical modes in the left and right cavity half (each 
of length respectively, whose resonance frequency 
ujo is changed due to the displacement x of the mem- 
brane. The coupling g describes the photon tunneling 
through the membrane. Due to the coupling, there is 
an avoided crossing in the optical resonance frequency 
oo±{x) = ±^/g^^^^J^JoxJ^)^ ^ see Fig. [Tj). We propose to 
drive the membrane with mechanical frequency Q and 
resulting amplitude A around a mean position xq. 



x{t) = Acos{Qt) + xo. 



(2) 



and investigate the system in the regime where the 
timescale of photon exchange is comparable to the 
timescale of the mechanical motion {g :^ Q). Recently 
the coupling frequency g/27r has been significantly re- 
duced by exploiting properties of transverse modes [20j, 
and it is tunable down to 200 kHz at present. The me- 
chanical eigenfrequencies of typical ImmxlmmxSOnm 
membranes range between 100 kHz and IMHz. Com- 
mercially available membrane sizes should allow this to 
go from 20 kHz up to 10 MHz. We point out that here 
Q need not coincide with the membrane's eigenfrequency 
but depends only on the driving. 

We assume the left hand side of the cavity is driven by 
a laser of frequency ujl and amplitude b^^.Om goal is to 
examine the photon dynamics by looking at the trans- 
mission T. Using input/output theory, the equations of 
motion for the average fields — {cll) and aR = {clr) 
read 



dt 
dt 



[-x{t) aL^gaR] - -aL - V^^r(t) 



(3) 



with the cavity decay rate n for each of the modes, and 
the drive ^^^(t) = ^-^^LtJJ^n^ Here, we used a rotating 
frame, with laser detuning from resonance Al = ool—^o- 
The displacement is written in terms of a frequency via 
x{t) = {u;o/l)x{t), likewise for A, xq. The transmission to 
the right, T{t) = n{a^^{t)aR{t))/{h'^)\ can be expressed 
as 



T{t) = 



(4) 



where the phase factor includes laser drive and cavity 
decay, while the Green's function G(t,t') describes the 
amplitude for a photon to enter the cavity from the left 
at time t' and to be found in the right cavity mode later 
at time t. Technically, G{t,t') is found by setting n = {) 
in Eq. (|3| and solving for aR{t) with the initial conditions 
aL{t') = 1, aR{f) = 0. 

We start investigating the dynamics by considering 
modest drive amplitudes A < Q. Fig. [2|l displays the 
time-averaged transmission depending on xq and A^. We 
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Figure 2: (a) Autler-Townes splitting of the cavity frequen- 
cies LJ± due to the mechanical motion: density plot for the 
time-averaged transmission depending on xq and A^. For 
every p osition x p the mechanical drive frequency is set to be 
Q = 2\/ + ^0 such that it is always resonant with the differ- 
ence between the two optical mode frequencies uj±. Further 
parameters are amplitude A — 0.2 Q and decay k — 0.1 5'. 
The splitting is set by the Rabi frequency g\ gA/Q, pro- 
portional to the drive amplitude, (b) Transmission as in (a) 
but for stronger drive A = 1.6 Mechanical sidebands, dis- 
placed by ±Q, become visible and interact. 



observe an Autler-Townes splitting [161 E] of the two 
hyperbola branches u;±. Indeed, the mechanical drive in- 
duces Rabi oscillations between the two photon branches, 
at a Rabi frequency gi c:^ gA/Vt, leading to a correspond- 
ing splitting in the spectroscopic picture. For larger drive 
amplitudes the dynamics becomes more involved. For in- 
stance mechanical sidebands arise as shown in Fig. [2]3, 
and they start to interact with each other. In the follow- 
ing, we will focus on the dynamics of this strong driving 
regime. 

Fig. [3] shows numerical results for A ^ ^, ^. For 
experimentally accessible parameters ^/27r = IMHz, 
/ = 1 cm and a;o/27r = 3-10^^ Hz, we have coq/I = 30^/nm 
and A = 60g corresponds to an oscillation amplitude of 
A — 2 nm that is below the nonlinear regime for a 50 nm 
thick membrane. Apart from the modulation of trans- 
mission as a function of A (see below), we observe finite 
transmission only if is a multiple of Q.. We first present 
an intuitive description. Transmission is determined by 
two subsequent processes. First, the laser has to excite 
the left mode. Secondly, the internal dynamics must be 
able to transfer photons into the right one. In general, 
both processes are inelastic and therefore require energy 
to be transferred between the light field and the oscillat- 
ing membrane. The left mode's frequency is oscillating 
around the time-averaged value — xq. Hence, the reso- 
nance condition to excite the left mode reads 



(5) 



see Fig. |4^. Here, m£t is an adequate multiphonon tran- 
sition. The width of the individual resonances is deter- 
mined by K. The subsequent process displays the physics 
of LZS dynamics: LZ transitions split the photon state 
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Figure 3: Density plot for the time- averaged transmission as 
function of average displacement xq and mechanical drive am- 
plitude A ^ g. Further parameters are laser detuning 
Al = 0, mechanical frequency Q/27r = 0.2g and cavity decay 
K — 0.2g. Finite transmission is observed when the resonance 
conditions Eq. ^ and ([6| for multiphonon transitions are 
met. The transmission is modulated according to the prod- 
uct of two Bessel functions. For the case xq = —Q both are 
depicted in the plot's plane. Red: ^ J^{A/Q), due to the ex- 
citation process. Yellow: ^ J^(2A/Q), due to LZS dynamics. 



into a coherent superposition, the two amplitudes gather 
different phases and interfere the next time the system 
transverses the avoided crossing. The condition for con- 
structive interference can also be phrased in terms of an 
additional multiphonon transition that transfers a pho- 
ton from the left mode with average frequency —xq to 
the right one at +xo, 



(6) 



We find transmission only if both conditions are met. We 
note that the coupling g between modes does not enter 
here. We will come back to this point later. 

To derive these resonance conditions as well as to un- 
derstand the dependence on A, in the following, we cal- 
culate an approximate, analytic expression for the trans- 
mission. From Eq. ([3|, the Green's function G(t,t'), re- 
quired for the transmission Q, is found to be 



(7) 



where we have split off a phase (l){t') = (A/Q) sm{Qt')^ 
and CLRit^t') is a solution to 
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with t>t' and initial condition dR{t' ^ t') = 0, dL{t\ t') = 
1. We now show that the two multiphonon processes in- 
troduced above correspond to the two factors in Eq. ^ . 
The term e-*^^^') = Y.^Jm{A/Vt)e-^'^^^' describes the 
initial excitation, where the amplitude for a transfer of m 
phonons is set by the Bessel function Jm{A/Q). Secondly, 
the internal dynamics described by dR{t^t') is expressed 
in terms of a two-level system with time-dependent cou- 
pling ^e^*^^^) = ^En -^n(2A/Q)e^^^^ Thus, the strength 
of the second multiphonon transition nVt in Fig. is de- 
termined by a Bessel function Jn{2A/Q)^ corresponding 
to Stueckelberg interferences known from atomic physics. 
As a special case, this also describes the Autler-Townes 
splitting at small drive. This can be calculated from ([8| 
using an interaction picture representation and consider- 
ing the time-dependent coupling only up to Ji, yield- 
ing an effective transition frequency 1\/ g\ + with 
go — gJo^AI^\ and a Rabi frequency g\ — gJxi^AjO). 

In the case of LZS dynamics, i.e. strong drive, for suf- 
ficiently large amplitudes only one of the harmonics of 
g^^Jn^Aj^^e^"^^^ will be in resonance with the sys- 
tem. This corresponds to leading-order perturbation the- 
ory within the Floquet approach [22 applied to Eq. ([8|. 
In this case Eq. ([s]) simplifies to the problem of a two- 
state system with harmonic drive at nVt and effective 
coupling constant 



gn=gJn{^Al^). 



(9) 



To estimate when this approximation becomes appropri- 
ate, we note that for a driven undamped two-state system 
the width of the power-broadened resonance is set by the 
Rabi frequency. Thus, Eq. ([8| yields a series of reso- 
nance peaks at xo = nQ./2^ and they become separated if 
Agn < Using the asymptotic form for large arguments 
A/^l > 1, J„(y) ^ y^cos (y - - f ), we find the 
resonance approximation to hold whenever 



^ 16 



(10) 



Note the resemblance to the criterion for non-adiabatic 
transitions that can be derived from the standard LZ for- 
mula Pi^i = exp(— 7r^^/2'y), where v = AVt is the sweep 



velocity. Eq. (10) is clearly fulfilled for the parameters of 
Fig. [31 

Given the resonance approximation, we find for the 
Green's function 

G{t,t') = -i^sin(c.;„(t-t'))e-*""(*+*')/2e-^*(*'), 

(11) 

with uon = \J {gnf + (^0 " nQ/2)'^. Note that cOn con- 
tains gn, which is much smaller than the bare splitting g 
for A ^ Q. This explains why the resonance conditions 
([5| and (|6| involve the bare optical mode frequencies ±xo 
instead of the adiabatic eigenfrequencies u;±. 
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We insert (11) into Q, taking into account the sum 
over independent contributions with n quanta. In the 
resolved sideband regime {Q > k)^ the integration of ffl 
selects a specific m for the excitation process, see Eq. (pF. 
Thus we find an approximate expression for the transmis- 
sion (displayed here for the special case Al = 0, where 
m = 2n): 



ALU - LUo 



(12) 



This captures fully the numerical results shown in Fig. |3] 
Whenever the resonance conditions are fulfilled (by 
choosing the offset Xq) , the transmission is modulated by 
the two Bessel functions. While Jm{A/^}) describing the 
excitation process depends on the amplitude the LZS 
dynamics characterized by Jni'^A/Q) is determined by 
the phase difference gathered between LZ transitions, in- 
volving 2 A. According to Eq. ([5|, if we were to increase 
Al in Fig. |3] we would tune out of resonance and the 
transmission would vanish everywhere. For Al = 1^/2, 
the conditions (|5| and ([6| can be met for xq being an 
odd multiple of Q/2. For Al = ft we would again find 
transmission for xq being a multiple of Q. Note however 
that the entire plot would be shifted in A by an amount 
7rft/2 due to the changed index of the Bessel function J^. 
Finally, at = 21] we would recover Fig. [3| 

Fig. [4JD shows numerical results for smaller values of 
A/g (while keeping Q as in Fig.jsj). As before, we see res- 
onances for xq being a multiple of ^} and expect regions of 
excitation with width determined by Jm{A/rt). Within 
these regions, we note the already familiar substructure 
that is due to LZS dynamics. 

To conclude, we proposed a setup to investigate non- 
equilibrium photon dynamics driven by mechanical mo- 
tion in an optomechanical system with a membrane in- 
side a cavity. We predicted the possibility to observe 
Autler-Townes splitting and features of Landau- Zener- 
Stueckelberg dynamics in the transmission spectrum. 
The observation of the effects discussed here is within 
reach of current experiments. The same nontrivial light 
field dynamics will enter when describing self-induced 
nonlinear optomechanical oscillations in these systems, 
which would be an interesting topic for future research. 

We acknowledge fruitful discussions with Jack Sankey, 
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